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Abstract 

In this paper we study the viability of persuing analytic variational techniques for the calculation 
of glueball masses in 3+1 dimensional Hamiltonian lattice gauge theory (LGT) in the pure gauge 
sector. We discuss the major problems presented by a move from 2+1 to 3+1 dimensions and 
develop analytic techniques to approximate the integrals appearing in 3+1 dimensional variational 
glueball mass calculations. We calculate ++ and l" 1 glueball masses on a lattice consisting 
of a single cube. Despite the use of a very simplistic model, promising signs of an approach to 
asymptotic scaling is displayed by the SU(iV) l" 1 glueball mass as N is increased. 
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I. OUTLINE 



In this paper we explore the viability of extending the analytic techniques used with 
success in Refs. Q and 3 to the calculation of glueball masses in the pure gauge sector in 
3+1 dimensions. The primary difficulty lies in the calculation of expectation values in 3+1 
dimensions. In Section|H]we briefly review what has been achieved in Hamiltonian LGT jjj in 
3+1 dimensions. We discuss the difficulties faced in 3+1 dimensions and possible solutions in 
Section UTTl In Section HVl we consider the problem of Gauss' law constraints. This is a topic 
that has been discussed in the context of Hamiltonian LGT most recently by Ligterink, Walet 

n 

and Bishop jj| and concerns the constraint equations that appear when non-abelian gauge 
theories are canonically quantised. In Section we move on to the calculation of variational 
glueball masses on a single cube using the analytic variational technique discussed in Ref. Q. 
We finish in Section IVll with a discussion of the viability of pursuing analytic techniques for 
pure SU(iV) LGT in 3+1 dimensions based on the results of Section IVl 

II. INTRODUCTION 

From a renormalisation point of view the key difference between 2+1 and 3+1 dimensional 
gauge theory lies in the units of the coupling constant. 2+1 dimensional gauge theory has 
a coupling constant, e 2 , with the dimensions of mass and so the coupling constant explicitly 
sets a mass scale for calculations on the lattice. In contrast the 3+1 dimensional coupling 
constant is dimensionless. This makes the extraction of continuum physics from lattice 
calculations more subtle in 3+1 dimensions than in 2+1. 

On a practical level, there is a more serious problem faced in moving from 2+1 to 3+1 
dimensions for Hamiltonian LGT calculations. The analytic techniques that were used with 
success m Refs. and 3 are no longer applicable. These techniques rely heavily upon the 
fact that in 2+1 dimensions a change of variables from links to plaquettes has unit Jacobian. 
The form for the equivalent Jacobian in 3+1 dimensions is considerably more complicated. 
The most comprehensive study of the change of variables from links to plaquettes is due to 
Batrouni P, |(| . We discuss this change of varia ales in more detail in Section IIII1 

The extension of the techniques used in Ref. yj to 3+1 dimensions is not straightforward. 
There are a number of immediate problems. Firstly, since plaquettes are not independent 
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variables in 3+1 dimensions one can not automatically work in the infinite volume limit. In 
a precise study one would need to calculate identical quantities on different sized lattices and 
extrapolate to the infinite volume limit. Secondly, in the context of analytic calculations, 
even on small lattices the integrals involved in the calculation of basic matrix elements 
are considerably more complicated than those encountered in 2+1 dimensions. Such matrix 
elements could in principle be carried out analytically on small lattices but since the number 
of integration variables increases quickly with the volume of the lattice a calculation on 
even a 5 3 lattice would seem exceedingly difficult. How quickly the infinite volume limit is 
reached will therefore determine the worth of pursuing analytic Hamiltonian methods in 3+1 
dimensions. Finally, there is the complication of Gauss' law which we discuss in Section II VI 

The only Hamiltonian techniques to have been applied with any success to the case of 
SU(3) gauge theory in 3+1 dimensions have been strong coupling expansions, the t-expansion 
and exponential wave function methods. Each of which we now summarise. 

Strong coupling perturbative techniques were used in the early days of LGT in an attempt 
to bridge the gap between the strong and weak coupling limits. Strong coupling expansions 
of the Callan-Symanzik (3 function were calculated to 0(g~ 2A ) [3, 0| and showed signs of 
interpolating smoothly between the strong and weak coupling limits. This suggested that 
the j3 function was a smooth function of the coupling with its only zero at g = 0, providing 
a strong argument at the time for the continuum limit of LGT confining quarks. Corre- 
sponding strong coupling expressions for glueball masses did not share the same success. 
Despite strong coupling calculations to 0(g~ 28 ) j^J scaling was not observed in ++ , l + ~ or 
2 ++ glueball masses. It was later realised that a roughening transition prevents a smooth 
crossover from strong to weak coupling physics (see Ref. Q and references within). 

The t-expansion was introduced by Horn and Weinstein as an analytic method suit- 
able for the study of LGT in the Hamiltonian formulation. It has been applied in the 
calculation of glueball masses in 3+1 dimensions for SU(2) [3| and SU(3) 14, 1^ LGT 
in the pure gauge sector. More recently it has been used in an attempt to calculate the 
lowest hadron masses [16]. For each case however asymptotic scaling of masses was not 
directly observed. Extrapolation techniques such as Pade approximants were required to 
probe the weak coupling region. The extrapolated mass ratio results agreed with Monte 
Carlo estimates of the time. 

The coupled cluster method and related exponential wave function techniques have re- 
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ceived by far the most attention in Hamiltonian LGT in recent years. Essentially these tech- 
niques aim to solve the Kogut-Susskind eigenvalue equation by making a suitable ansatz for 
the wave function. The coupled cluster method was originally constructed with applications 
in nuclear physics in mind [l?! Q] but has since found the majority of its applications in 
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molecular ph ysic s |ig| | . Its application in the context of Hamiltonian LGT is described in 



Refs. |20| and |21|. The truncated eigenvalue method, developed by Guo, Chen and Li 
is another exponential wave function technique to have found application in Hamiltonian 
LGT. 

A number of groups have made considerable progress in the application of exponential 
wave function techniques to Hamiltonian LGT in 3+1 dimensions. While most studies have 
explored gauge groups other than SU(3) in less than three dimensions, studies of SU(3) 
glueballs in 3+1 dimensions have commenced. Results tangent to the expected scaling 
form, indicating an approach to scaling, have been obtained for SU(3) pure gauge theory in 
a coupled cluster calculation by Leonard [23]. However convergence with increasing orders 
appears to be slow. Finite order truncation errors appear to be under more control in 
the truncated eigenvalue method. The first calculations of 3+1 dimensional SU(3) glueball 

n 

masses with this method |24| gave a 1 + to ++ mass ratio that was consistent with the 
Monte Carlo results of the time. The agreement was not as good for the to ++ mass 
ratio. A convincing demonstration of asymptotic scaling has not yet been produced in a 
calculation of glueball masses for SU(3) in 3+1 dimensions. 

More promising 3+1 dimensional results have been obtained for higher dimensional gauge 
groups. Ironically these results have been obtained with much simpler methods than either 
of the exponential wave function methods described above. Asymptotic scaling of the lowest 
++ glueball mass has been demonstrated by Chin, Long and Robson in a variational calcu- 
lation on a small volume (6 3 sites) lattice for SU(5) and SU(6) This calculation follows 
the variational technique described in Ref. Q but uses Monte Carlo rather than analytic 
techniques to calculate the required expectation values. It uses only plaquette states in the 
minimisation basis rather than the large basis of rectangular states used in Refs. QandQ 
Naively one would hope that the same method could be explored on similar sized lattices 
for larger N, with additional states in the minimisation basis, using the analytic techniques 
of Ref. Q. We take this as our motivation for the studies presented in this paper. It will 
become clear however that attempting a similar calculation to Chin, Long and Robson on a 
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single cube using analytic techniques presents a significant challenge. 



III. THE MOVE TO 3+1 DIMENSIONS 

In LGT one usually works in 2+1 dimensions to test a technique with the intention of 
later extending it to the physically relevant 3+1 dimensions. This may be justified in the 
Lagrangian approach where the time coordinate is treated on the same footing as the spatial 
coordinates and adding another dimension equates to nothing more than an increased load 
on computer memory. However, in Hamiltonian LGT two serious technical differences exists 
between 2+1 and 3+1 dimensions. The first is Gauss' law and the second is related to 
constructing a Jacobian for transforming from link to plaquette variables. 

To understand the first difference one needs to recall that Hamiltonian LGT is formulated 
in the temporal gauge which sets A = 0. If one starts with the Yang-Mills Lagrangian and 
performs the standard equal-time quantisation, one runs into problems because the time 
derivative of Aq does not appear in the Lagrangian. The variational principle gives equations 
of motion for the space-like components resulting in the standard Yang-Mills Hamiltonian. 
For the time-like component we obtain a set of algebraic constraint equations which are the 
analogue of Gauss' law. There is one constraint equation for each colour component of A , 
N 2 — 1 for SU (N), at each lattice site. As has been pointed out in the context of Hamiltonian 
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LGT by Ligterink, Walet and Bishop |4J, it is only when one works with a set of variables 
whose number of degrees of freedom matches the number of unconstrained degrees of freedom 
in the theory, that one can avoid the technicalities of constraint equations. For this reason, 
in 2+1 dimensions we really are quite lucky. The number of plaquette variables on a square 
two dimensional lattice is precisely equal to the number of unconstrained variables [4|. For 
a three dimensional cubic lattice this is not the case. One could envisage constructing a 
polyhedral lattice in d dimensions such that the number of faces (plaquettes) would equal 
the number of unconstrained variables. However such lattices appear to be prohibited by 
Euler's equation which relates the number of vertices, edges and faces of polyhedra. 

The other serious technical difference between 2+1 and 3+1 dimensions, that of con- 
structing a Jacobian for transforming from link to plaquette variables, is relevant to both 
the Hamiltonian and Lagrangian formulations of LGT. It only becomes important when the 
method of choice relies on plaquette variables for its implementation. The transformation 
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from link variables to plaquette variables is required in the approach taken in Refs. 
to make use of the analytic results available for certain group integrals. Such a transforma- 
tion does not need to be made if one is happy to use Monte Carlo techniques to handle the 



integrals as is the case in Ref. 
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The most complete treatment of the transformation from link to plaquette variables is 
due to Batroun, fifl who worked in the L„» foliation. His approach was based 
on the continuum work of Halpern 26] who constructed a field-strength formulation of gauge 
theory. In Halpern's construction the definition of the field strength is inverted to give 
an expression for as a function of the field strength. To do this requires the choice of 
a suitable gauge. The Jacobian of the transformation is precisely the Bianchi identity; a 
constraint equation on F^ u . For the Abelian case the Bianchi identity is equivalent to the 
requirement that the total magnetic flux leaving a volume is zero. Batrouni developed an 
equivalent field strength formulation on the lattice. On the lattice, the link elements, C/j, 
correspond to the vector potentials, A^, and the plaquette variables correspond to the field 
strengths, F^. Batrouni demonstrated that the Jacobian of the transformation from link 
variables to plaquette variables was the lattice analogue of the Bianchi identity. For Abelian 
gauge theories the lattice Bianchi identity can be separated into factors which depend only 
on the plaquette variables of elementary cubes, with one factor for each cube of the lattice. 
For non-abelian theories the Bianchi identity has only been found to separate in this way 
for special types of lattices, the largest volume example being an infinite tower of cubes. 
For the infinite lattice the Bianchi identity is a complicated nonlinear inseparable function 
of distant plaquette variables. Interestingly it is the only source of correlations between 
plaquette variables in LGT. Mean plaquette methods have been developed to deal with the 
added complications of the non-abelian Bianchi identity but have not progressed far 27 1. 

To summarise, Hamiltonian LGT in 3+1 dimensions faces some serious problems. Care 
must be taken in choosing appropriate variables to work with if constraints on the lattice 
electric fields are to be avoided. Additionally, if one wishes to use analytic techniques to 
calculate matrix elements the complications of the Bianchi identity restricts the calculations 
to small lattices. 
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IV. CONSTRAINT EQUATIONS 



In this section we focus on one of the problems faced in moving from 2+1 to 3+1 di- 
mensions, that of Gauss' law. From the discussion of Section 11111 it would seem that one 
needs either to construct an appropriate set of variables with precisely the correct number 
of unconstrained degrees of freedom or be faced with the problem of building constraint 
equations on the lattice electric fields. As has been pointed out by Ligterink, Walet and 
Bishop jj] there is an alternative. The problem of satisfying Gauss' law can be solved by 
working with wave functions that are annihilated by the generator of Gauss' law. We discuss 
this matter in what follows. 

The generator of Gauss' law on the lattice can be written as P, Q , 

g*(x) = ^2[£r(x)+S^(x)]. (1) 

i 

Here Sf(x) is the lattice chromoelectric field on the directed link running from x to x + ai. 
In this notation the lattice electric fields satisfy the commutation relations l| 

[£^x),Uj(y)} = T*U j (y)5 1d 5 ay (2) 
Sliix + ia^Ujiy)] = -U^T^jSny. (3) 

Here i is a unit vector in the i direction and {T a : 1 < a < N 2 — 1} is a basis for SU(iV). It is 
common to use the Gell-Mann basis, in which case T a = A a /2, where {A a : 1 < a < N 2 — 1} is 
the set of traceless N x N Gell-Mann matrices. It should be pointed out that in this notation 
we have U^(x) = U^i(x + ia). For physical states we must therefore have Q a (x)\ip) = for 
each lattice site, x, and all a = 1, ... , iV 2 — 1. It should be checked that this is the case for 
the one plaquette exponential trial state, 

|0 o >=exp(c^ReQJ |0). (4) 

\ x,i<j ) 

Here, |0) is the strong coupling vacuum defined by S°"(x)\0) = for all i, x and a = 
1,2,..., iV 2 — 1. The directed square denotes the traced ordered product of link operators, 
Ui(x), around an elementary square, or plaquette, of the lattice, 

Q := Tr [Ui(x)Uj(x + ia)U}{x + ja)U}(x)] ( 5 ) 
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where a is the lattice spacing. 

Consider first a state consisting of a single plaquette acting on the strong coupling vac- 
uum, 

\pij(x)) = Tr [[/;(:*;) £/i(cc + ia)U^(x + ja)U^(x)j |0>. (6) 
Since |0) is annihilated by the electric field, using Eqs. (j2J) and © we immediately have 



g a {v)\p td {x)) = o, 



(7) 



for all lattice sites, y, not lying on the corners of the plaquette, Pij{x). Consider now sites 
lying on the corners of the plaquette in question. In particular consider y = x. Making use 
of the commutation relations of Eqs. (J2J) and (jHJ) we have 

g a (x)\ Pij (x)) = Tr\T a U i (x)U j (x + ta)U^(x + ja)U-. j (xj\ |0) 



-Tr 



Ui{x)Uj(x + m)U^(x + ja)U^(x)T a |0) 



0. 



(8) 



The same applies for other sites on the plaquette. This argument is not specific to plaquettes. 
Gauss' law is found to be satisfied locally by any closed Wilson loop on the lattice, traced 
over colour indices, acting on the strong coupling vacuum. It is easy to extend this result 
to products of such loops acting on |0). To see this we consider how Gauss' law applies to 
the product of two plaquettes 



\Pij(x)pij(x + at)) = Pij(x) Pij (x + ai)\0). 



(9) 



Only at the sites x + ai and x + a(i + j) does this case differ from the single plaquette 
example . Let us consider Gauss' law at x + = x + ai. Once again, using the fact that the 
strong coupling vacuum is annihilated by the lattice chromoelectric field, and Eqs. Q and 
(jSJ) we have 

g a (x + )\ PlJ (x) PlJ (x + )) = ^2[S^{x + )+£l k {x + ), Pij {x)p ij {x + )] |0> 

k 

= { [ £a kM + £t k (x+), Pij (x)] Pij (x + )\o) 

k 

+ Pij (x) [S%{x + )+Sl k (x + ), Pij {x + )] |0>} 
= 0. (10) 



This result is easily extended to arbitrary products of closed Wilson loops, traced over colour 
indices, acting on |0). Consequently, any function of such loops acting on |0) also satisfies 
Gauss' law provided it admits a Taylor series expansion. The one plaquette trial state of 
Eq. (@J) thus obeys Gauss' law and is therefore suitable for use in simulating physical states. 



V. THE ONE CUBE UNIVERSE 

In this section we adapt the analytic techniques used with success in Refs.QandS for use 
in 3+1 dimensions. As a starting point we consider the case of a single cube. This will serve 
as a test case to assess the feasibility of a larger volume study. Our aim is to use a small 
basis of states to calculate 3+1 dimensional SU(iV) glueball masses variationally, and check 
if an approach to the correct scaling form is observed. We do not expect to achieve scaling 
with such a small lattice. However, an approach to scaling would warrant an extended study 
on larger lattices. We start with a general description of the approach. The starting point 
is the choice of trial state. As in 2+1 dimensions we choose to work with the one 

plaquette trial state of Eq. (jlj). 

Many Hamiltonians are available for LGT calculations j^ . In Ref. Q we used Kogut- 
Susskind |3j, improved and tadpole improved Hamiltonians to calculate glueball masses for 
N < 5. Here we employ the simplest of these, the Kogut-Susskind Hamiltonian, which is 
defined for pure SU(iV) gauge theory with coupling, j 2 , on a lattice with spacing, a, by, 

x,i fi*<i 

where the plaquette operator is given by, 

Piii*) = 1-^ReQ. (12) 
The variational parameter, c, is fixed by minimising the vacuum energy density, 

« = ±m 

N 2 -l\ 201 /r-i\ _„ ( 13 ) 



2/3 
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as was done for the case of 2+1 dimensions in Ref. fl. Here (3 = N/g 2 and N p is the number 
of plaquettes on the lattice. In 3+1 dimensions, however, the expression for the plaquette 



ueball mass 
l[ The key 



expectation value in terms of the variational parameter is significantly more complicated than 
for the case of 2+1 dimensions. Having fixed the variational parameter we then construct a 
small basis of states, with each state fitting on a single cube, and minimise the g 
over this basis. The process of minimising the massgap follows precisely Ref. 
difficulty of working in 3+1 dimensions is the calculation of the required integrals. In the 
next section we explain how an analytic approximation to these integrals can be obtained. 

A. SU(iV) integrals in 3+1 dimensions 

The analytic techniques used in Ref. [l rely on the fact that the transformation from link 
to plaquette variables has unit Jacobian in 2+1 dimensions. Batrouni has calculated the 
Jacobian for arbitrary numbers of dimensions 0,0]. For a lattice consisting of a single cube, 
the result of Ref. y is 

J = 5{P 1 P 2 P i P^P & - 1) 

r r i=l 

where Pi, . . . ,P e are the six plaquette variables on the single cube and the sum is over all 
characters, Xr, of SU(iV). d r denotes the dimension of the character, \ r . The second line is 
simply a character expansion of the first line. In the variational study of glueball masses, we 
need to calculate the integrals of overlapping trace variables on a single cube. It is always 
possible to reduce these integrals to integrals involving non-overlapping trace variables using 
the orthogonality properties of the characters. For example, consider the expectation value, 
on a single cube, of a twice covered bent rectangle, with respect to the one-plaquette trial 
state of Eq. (@J), 

Y r ^ « (15) 



To proceed with this integral we need to perform character expansions. The orthogonality 
of the characters can then be used to write characters over two plaquettes in terms of one 
plaquette characters. To demonstrate how this is done we consider the character expansion 
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of one of the integrals in Eq. (|15|) , 

J dP Xr {P) X2 {PP>)e^ p+p ^ = Y^c r , J dP Xrl {P)X2{PP') 
Here the <v are given by 

c r , = J dP X AP)Xr(P)e^ P+P ^. 
Making use of the orthogonality property of characters given by, 

J dU p Xr'{U p V)Xr{U p ) = ^-5r'rXr(V), 



(16) 



(17) 



(18) 



we obtain 



dP Xr (P) X2 (PP')e cTr(P+Pt) = ^ X2 (P) X2 (PP') 

d 2 



(19) 



Making use of this, and proceeding similarly for the analogous integral involving XuiPsPe 
we can reduce Eq. (|T5j) to 



d 2 



^ |^o> = E^[/ dPXr{P)e™^ 
/ dQ X 2(Q)Xr(Q)e^ Q+Q ^ + dQ Xll (Q)Xr(Q)e^ Q+Q ^ 



(20) 



The plaquette matrix element provides a more straightforward example; 

(^o|Q|0o> = £jf [/ dP X r{P)e^ P ^ J dQTiQxr(Q)e cTr{Q+Q ^ 



(21) 



To proceed further we need expressions for the integrals over characters that appear in each 
matrix element of interest. We define these character integrals generically by 



C ri r2—r„ ( c ) 



J dP Xri (P)Xr 2 (P) ■ ■ ■ X^)e cTr(P+Pt) . 



(22) 



All SU(iV) integrals encountered in the calculation of glueball masses in 3+1 dimensions 
can be expressed in terms of character integrals. It is possible to calculate them generally, 
however three points need to be considered. Firstly, if we are to use the machinery of Ref. Ill, 
we need to know how to express a given character in terms of trace variables. Secondly, 
as the dimension of the gauge group increases the number of independent trace variables 
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increases rapidly. In order to contain the number of integrals required for a given order of 
approximation, it is necessary to express all high order trace variables in terms of a basis 
of lower order ones. Finally, the collection of SU(iV) integrals presented in Ref. must be 
extended. We now address each of these points in turn. 

The problem of expressing the characters of SU(iV) in terms of trace variables was solved 
by Bars 29] who showed that for SU(iV), 

Xn(U) = £ 5 (X> -n) f[ * (TH/*)* (23) 

ki,...,k„ \i=l / j=l j ' J 

General characters can then be expressed in terms of Xn{U) using standard techniques from 
group theory j^], 

Xnr 2 -r N - 1 (U) = det [xr J+ i-j(U)] 1 < ij < iv _ 1 • (24) 

Here the quantities inside the determinant are to be interpreted as the (i,j)-th entry of 
an N x N matrix. For a given character, Eq. (|2*4*j) produces a multinomial of traces of 
different powers of U. However, not all of these trace variables are independent. The so 
called Mandelstam constraints define the relationship between dependent trace variables. 
For SU(2) and SU(3) we have 

TrU = TrC/t W G SU(2) and 

(25) 

Tr(U 2 ) = (TrU) 2 - 2Tr(C/t) WU G SU(3). 

Similar relations can be obtained to express all trace variables in terms of TrU for SU(2), 
and for SU(3) in terms of TrU and TrW . To reduce the calculation of character integrals 
to a manageable size, we require expressions for high power trace variables in terms of a 
minimal set of lower order trace variables. To do this we proceed as follows. We start with 
an alternative expression of det U = 1, satisfied for all SU(iV) matrices U, 

Here the colour indices of the group elements have been made explicit and all repeated 
indices are summed over, e^...^ is the totally antisymmetric Levi-Civita symbol defined to 
be 1 if . . . , i n } is an even permutation of {1,2, ... ,n}, —1 if it is an odd permutation 
and otherwise (i.e. if an index is repeated). Multiplying Levi-Civita symbols produces 
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sums of products of delta functions, the precise form of which depends on how many pairs 
of indices are contracted. A standard result from differential geometry will be useful here, 

e ai ...a n e bl ... bn = det (X^ )^ ^ • (27) 
Multiplying both sides of Eq. (ffilj) by Uj N j N+i and contracting over repeated indices gives 

e iii2—iN-ljN+iUi 1 j 1 Ui 2 j2 ' ' ' Ui N _ 1 j N _ 1 = £ j 1 j2—jN^j N j N+1 
e hh---lN-iiN e hi2---iNUi-ij-iUi 2 j 2 ■ ■ ■ Ui N _ 1 j N _ 1 = £l 1 l 2 ...l N ^ii N £ jij2--jN^j N i N - (28) 

In the last line we have renamed dummy indices and multiplied both sides by a Levi-Civita 
symbol. Making use of Eq. (J27)l in Eq. (J2HJ) produces a product of delta functions on each side 
of the equation. By introducing appropriate combinations of delta functions and contracting 
over repeated indices we can construct trace identities as we please. We find that identities 
useful in our context are produced by the introduction of delta functions and an additional 
matrix, A, which is not necessarily an element of SU(iV), into Eq. (|28|) as follows: 

N-l N-l 
£k 1 —k N - 1 i N £i 1 —i N Ui 1 j 1 Aj 1 k 1 J^J fik m j m Ui m j m = £k 1 —k N - 1 i N £j 1 —j N Aj 1 k 1 Uj NiN J^J fikmjm- (29) 

m=2 m=2 

Here we have again renamed dummy indices. Some examples of SU(iV) identities obtained 
from Eq. (J29j) valid for all N x N matrices A are as follows: 

- TcUTc(UA) + Tr(U 2 A) = -TrATrf/ f + Tr(U^A) W G SU(3) 
-(TtU) 2 Tt(UA) + Tt(UA)Tt(U 2 ) + 2TrUTi(U 2 A) - 2Ti(U 3 A) 

= -2TrATrC/ t + 2Tr([/ t A) W G SU(4). (30) 

We require formulae for Tr(U n ) in terms of lower order trace variables for SU(iV). Such 
formulae, known commonly as Mandelstam constraints, are obtained by setting A = Jj n - N + 1 
in Eq. ([29)1 . In Appendix El we list the Mandelstam constraints obtained in this way up 
to SU(8). We see that for SU(iV) it is possible to express all characters in terms of N — 1 
trace variables. In what follows we will choose to express the general SU(iV) characters in 
terms of the set of trace variables {TrW, TrU, TrU 2 , . . . , TyU n ~ 2 }. It may well be possible 
to reduce the size of this set. This would indeed improve the efficiency of our technique. 
However, no effort has been made to do this at this stage. 

The procedure for calculating the general character integral C ri ___ rn (c) is then as follows. 
We first express the characters Xn , ■ ■ ■ > Xr„ in terms of trace variables using Eqs. (|2*3*|l and 
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). We then simplify these expressions using the Mandelstam constraints of Appendix 1X1 
The final step is to perform the integrals over trace variables, which is an increasingly 
non-trivial task as one increases the dimension of the gauge group. For instance, the 
general character integral for SU(iV) involves integrals over powers of the trace variables 
TrW , TrU, TrU 2 , . . . TrU N ~ 2 . To proceed we need the following integral, 

T^- 9 f(c)= / dU(TiU qi ) Pl (TrU q2 ) P2 ■■■(TiU qk ) Pk e cTv{u+ut) . (31) 

JsV(N) 

The cases of interest to us here for SU(iV) are described by, T^_Yiz^7n-2( c )- To calculate 
this integral we need to extend the work of Ref. Q and consider the generating function, 



G 9 ,.. gfc (c, 7l , . . . , 7fc ) = / dUe ^u + u^U^^)_ (32) 
Following the procedure of Ref. Q we obtain, 

oo 

G qi ... qk (c, 71, . . . , Tfc) = ^ det [\i+j-i iqi: ... tqk (c, 71, . . . , 7fc)]i<y<jv > 

(33) 



with, 



A m , ?1 ,...,g fe (c,7i, ...,7k) = 2^ - I m+Siqi+ ... +Skqk (2c). (34) 

51,53,. ..,5^=0 

We then obtain an expression for T q Pl ; q p k k (c) by differentiating G gi _ qk appropriately; 

k 

(35) 



Qpi-\ hPfe 

^qi...q k k ( c ) = Q^px ; ; ; G gi ... gk (c, 71, . . . , 7^) 



7l=-=7fc=0 



B. The variational ground state 

In this section we fix the variational ground state following the usual procedure of min- 
imising the unimproved vacuum energy density given by Eq. (|13j). This equation is inde- 
pendent of the number of dimensions. The dimensionality of the lattice arises at the stage 
of calculating the plaquette expectation value. 

Making use of Eq. ([21)1 . the variational parameter can be fixed as a function of f3 for the 
one-cube lattice. In practice, the character sum in Eq. (|21|) and the infinite /-sum in Eq. (I33|) 
need to be truncated. We truncate the infinite /-sum at ±l max and instead of summing 
over all SU(iV) characters, we sum over only those characters, r = (r 1 ,r 2 , . . . ,r N _i), with 
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FIG. 1: The one cube SU(iV) variational parameter as a function of (3 for various N showing the 
dependence on the character sum truncation. 

r max > r i > r 2 > ' ' • > r iV-i- With this truncation scheme, memory constraints restrict 
calculations of the variational SU(iV) ground state to N < 7, when working with r max = 2 
and Imax = 2 on a desktop computer. 

The dependence of the variational parameter on r max with l max = 2 is shown for various 
gauge groups in Fig. [TJ As r max increases the variational parameter appears to converge 
for each N considered. Moreover the convergence appears to improve as the dimension of 
the gauge group is increased. With the exception of SU(3), the r max = 1 and 2 results are 
indistinguishable on the range < f < 0.7, where f = l/(Ng 2 ). As this is the ranee of 
interest to us later in the paper we restrict further calculations to r max = 1. 

The SU(iV) variational parameters on the one-cube lattice with r max = 1 and l max = 2 
are shown for various iV in Fig. El The results do not differ greatly. The corresponding 
variational energy densities are shown in Fig. El 
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FIG. 2: The one cube SU(iV) variational parameter in units of N as a function of l/(Ng 2 ) for 
various N. 
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C. Expressions for the glueball mass 

We follow precisely the method described in Ref. for the calculation of 3+1 dimensional 
glueball masses on a single cube. Here, however, we choose a different basis of states to 
minimise over. Instead of rectangular loops, we use states which fit in a single cube. We 
start with a basis of two states, the plaquettes and the bent rectangles, 



B = {|1>,|2», 



(36) 



with 



(37) 



and 



F 1 (x) = 
F 2 (x) = 



± 



(38) 



In 3+1 dimensions, the "+" sign corresponds to the ++ state and the "— " sign corresponds 



to the 1 + state 0. In order to calculate the glueball masses, following Arisue [31], we need 



expressions for the matrix elements N$, and D^ v 

1 



131], 



N r , 



Jj defined by, 
(l\H-E \l') 



(39) 



and 



(l\l f ) = J2 \(F}(x)F v (0)) - (FiWiFviO)) 



(40) 



Here the superscript, C, denotes the charge conjugation eigenvalue, C = ±1, of the state in 



question. Following Arisue 



improved Hamiltonians in Ref. |l|) : 



1 



we have the following reduction (which was generalised to 



(41) 
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After carefully counting the number of possible overlaps between different loops we arrive 
at: 




Here we have introduced the notation 

«0» = (0o|O|0o>. (43) 

The combinatorics which lead to the coefficients in the matrix elements in Eq. (142)1 are a 
result of counting the possible overlaps within a single cube. The remaining matrix elements 
N 2 2 and _D 2 2 can be calculated similarly, resulting in 




Having calculated N^, and D?,, we follow the minimisation procedure described in Ref. [l|to 
arrive at the glueball mass, AM PC . 
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D. Results 



In this section we present calculations of the ++ (symmetric) and l + ~ (antisymmetric) 
glueball masses on the one-cube lattice for SU(iV) with 4 < N < 7. We first define the 
rescaled glueball masse, fi PC , corresponding to AM PC as follows: 

f° = log(aAM^r 51/m ) - ^ log (^f) + (45) 

Here, as usual, P and C denote the parity and charge conjugation eigenvalues of the state 
in question. As a result of standard renormalisation group arguments 0, 3^|, asymptotic 
scaling of a glueball mass is observed if the corresponding rescaled glueball mass becomes 
constant for some range of couplings. 

The results for the rescaled symmetric glueball mass are shown in Fig. 0] Fig. |4(a)| 
shows the rescaled glueball mass calculated with only plaquettes in the minimisation basis. 
Fig. |4(b)"| shows the same quantity calculated with the minimisation basis of Eq. (|36|) . The 
aim of this exploratory study is to observe whether or not a move towards scaling is apparent 
as the number of states in the minimisation basis is increased. This is clearly not the case 
for the symmetric glueball mass on the range of couplings explored here. There is no sign 
of n ++ becoming constant on any range of couplings within < £ < 0.55. 

The prospects are marginally better for the antisymmetric case. The results for the 
rescaled antisymmetric glueball mass are displayed in Fig. A move towards scaling is 
observed, most clearly for the N = 7 case, near £ = 0.38. This matches the region of 



couplings for which Chin, Long and Robson observed scaling of the ++ glue 
N = 5 and 6 on a 6 3 lattice using only plaquettes in their minimisation basis 



>a 



1 mass for 



E. String tension 

If the exploratory study presented here was to be extended to larger lattices, at some 
stage it would be useful to compute the string tension, a. It is common for calculations in 
Lagrangian LGT to express results for masses in units of y/a. This allows masses calculated 
on the lattice to be expressed in MeV since the string tension can be calculated in MeV 
from the decay of heavy quarkonia for example. In the Hamiltonian formulation, precise 



calculations of the string tension have on 
Variational estimates, at least for SU(2) 



y 



Deen performed in the strong coupling regime. 
, have not exhibited asymptotic scaling when 
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FIG. 4: The one cube SU(iV) ++ rescaled glueball mass as a function of l/(Ng 2 ) for various N 
obtained with different minimisation bases. 

making use of the one plaquette ground state of Eq. (jlj). It is possible that the situation 
may be improved for higher dimensional gauge groups but this has not yet been tested. 

In this section we calculate the symmetric and antisymmetric glueball masses in units of 
y/a. For the string tension, since reliable variational results are not available, we use the 
strong coupling expansions of Kogut and Shigemitsu 34]. These are available for SU(5) and 
SU(6), as well as SU(2) and SU(3) which we do not consider here. 

Since the square root of the string tension has units of mass, the ratio of a mass to the 
string tension is constant in a scaling region. Again, with the crude model presented here, 
we do not expect to observe scaling. We seek only an indication of an approach to scaling. 
Such behaviour would warrant further study. 

The results for the symmetric states are shown in Fig. |U1 Since we perform calcula- 
tions with two states in the minimisation basis, the two lowest mass states are accessible. 
In Fig. |6(a)| the results for the lowest mass state are shown. Calculations with different 
minimisation bases are shown, with the "SU(iV)-plaquettes" label indicating that only pla- 
quette states are used. For each case the glueball mass has a local minimum in the range 
0.33 < £ < 0.37. The various results do not differ greatly. There is no sign of improved 
scaling behaviour when either bent rectangles are included or N is increased. The effect of 
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FIG. 5: The one cube SU(iV) 1 + rescaled glueball mass as a function of l/(Ng 2 ) for various N 
obtained with different minimisation bases. 

including bent rectangles is to lower the local minimum. 

The second lowest glueball mass is shown in Fig. 6(b) Again, no improvement in scaling 



behaviour is seen as N is increased. In Fie. |H1 the horizontal lines indicate the 3+1 di- 
— SU(5) calculations (and error bars) of Lucini and Teper Q. Interesting., both 
masses calculated here with a simplistic model, have minima that lie within the error bars. 

We now move on to the antisymmetric states. The results for the antisymmtric glueball 
masses in units of \fo are shown in Fig.0 The results are more promising than the symmetric 
case. We see an improved approach to scaling when bent rectangles are included in the 
minimisation basis for the lowest glueball mass. For this case the glueball mass shows 
promising signs of becoming constant in the ranges 0.32 < £ < 0.39 for SU(6) and 0.3 < £ < 
0.37 for SU(5). Only marginal improvement in scaling behaviour is evident as N is increased 
from 5 to 6. For each of these cases no data is available for comparison to our knowledge. 
Such promising signs are not apparent in the second lowest glueball mass in this sector as 
seen in Fig. |7(b)[ 

It would be interesting to perform analogous calculations for SU(7), for which the most 
promising results were displayed in Section IV Dl However strong coupling expansions for 



the SU(7) string tension are not available in Ref. 



34 



and have not been published elsewhere 
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FIG. 6: The one cube SU(5) and SU(6) ++ glueball masses in units of y/a as functions of 1/ (Ng 2 ) 
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to our knowledge. 
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VI. FUTURE WORK 



In this paper we have studied the variational ++ and l +_ glueball masses on a single 
cube in 3+1 dimensions. Our intention was to determine the viability of using analytic 
Hamiltonian methods in 3+1 dimensional glueball mass calculations. With such a crude 
small volume approximation and a minimisation basis containing only two states, the obser- 
vation of asymptotic scaling was not expected. Promising signs of an approach to asymptotic 
scaling was displayed by the l + ~ glueball mass, as N was increased, at couplings close to 
the scaling window observed by a comparable, but larger volume, study by Chin, Long and 
Robson [25]. Interesting results were also observed for calculations of the glueball masses 
in units of the string tension. With no variational results for the string tension available, 
the strong coupling results of Kogut and Shigemitsu [j^ were used. Interestingly, the mass 
of both SU(5) ++ glueball states calculated had minima which were consistent with the 

n 

Lagrangian calculation of Lucini and Teper [35J. Better scaling behaviour was exhibited 
by the SU(5) and SU(6) l + ~ states, although no alternative calculations are available for 
comparison to our knowledge. For the lowest mass l + ~ states calculated, the scaling be- 
haviour was improved by increasing the number of states in the minimisation basis and the 
dimension of the gauge group. The promising results observed in this paper warrant further 
study on larger lattices, with additional states in the minimisation basis. 

To extend this calculation to larger lattices and minimisation bases, two challenges will 
be faced. Firstly, the correct implementation of the Bianchi identity, within our plaquette 
based analytic approach, forces the number of integration variables to grow quickly with 
the volume. New approaches to handling the general character integrals may need to be 
developed in order to avoid the inevitable memory restrictions. Secondly, great care will 
need to be taken in the counting of overlaps between different states. This process could, in 
principle, be automated using techniques from graph theory and symbolic programming. 

It would also be of interest to extend the calculation presented here to larger N. To do 
this a more efficient technique for handling the general character integrals would need to 
be developed to minimise the memory demands of the calculation. The following is just 
one possibility. It is likely, that by expressing the products of characters appearing in the 
character integrals as sums of characters, one could reduce the memory requirements of the 
calculation drastically. The general character integrals could then be expressed in terms of 
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integrals of the form, 

G r (c) = [ dU Xr {U^ u+ur \ (46) 

JSXJ(N) 

where r = (r l5 r 2 , . . . r^r-i) labels a representation of SU(iV). Having completed the cal- 
culations presented in this paper, it was discovered that G>(c) can be handled using the 
techniques of Ref. Q, with the result, 

oo 

(47) 

l=— oo 

The final stage of this improvement would be to find a convenient way to express the products 
of characters appearing in the calculation as a sum of characters. This would be possible 
with the symbolic manipulation of Young tableaux. 
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APPENDIX A: MANDELSTAM CONSTRAINTS 



TrU n = TrU n - 3 - TrWTrU n - 2 + Trf/Trf/™" 1 W E SU(3) (Al) 



TiU n = -TiU n ~ 4 + Tr^TrC/"- 3 - ^(TiU) 2 TrU n - 2 + ^TiU 2 TrU n ~ 2 

+TrUTrU n - 1 VU E SU(4) (A2) 



TrU n = TrU n - 5 - Tr^TrfT" 4 + ^(TiU) 3 TiU n ' 3 - ^Tr^/Tr^W 1 " 3 

3 2 K ' 2 

+TrUTrU n - 1 VU E SU(5) (A3) 
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TrU n = -TrU n - 6 + TrU^TrU 71 ' 5 - — (Trin 4 Tr[/ n - 4 + -(TyU) 2 TiU 2 TrU n - i 

24 4 

-i(Trf/ 2 ) 2 Trf/ n " 4 - ^TiUTiU 3 TTU n ' 4 + ^TyU^yU' 11 ' 4 + i(Trf/) 3 Trf/ n " 3 

--TiUTrU 2 TrU n - 3 + -Tr£/ 3 Tr£/ n ~~ 3 - -Tr[/ 2 Tr£r~ 2 + -TrU 2 TrU n ~ 2 
2 3 2 2 

+TrC/TrC/ n - 1 VC/ G SU(6) (A4) 

TrC/ n = TrU n - 7 - TrC^lYC/ 71-6 + — (TrC/) 5 TrC/ n_5 - — (TrC/) 3 TrC/ 2 TrC/ n - 5 

120 v ; 12 v ; 

+ iTrf/(Trf/ 2 ) 2 Trf/ n - 5 + i(Trf/) 2 W 3 W n - 5 - ^TrC/ 2 TrC/ 3 TrC/ n ~ 5 

-iTrf/Trf/ 4 Trf/"" 5 + -TrU 5 TrU n - 5 - — (TrC/) 4 TrC/ n " 4 + i(Trf/) 2 Trf/ 2 Tr[/ n - 4 
4 5 24 4 

-i(W 2 ) 2 Trf/ n ^ 4 - iTrf/Tr^Trf/™- 4 + ^TrU 4 TrU n ~ 4 + jj;(TrC/) 3 TrC/ n - 3 

-iTrf/Trf/ 2 Trf/"~ 3 + -TyU^TyU^ 3 - ^-(TrC/) 2 lYC/ n_2 + ^TrC/ 2 TrC/ n_2 
2 3 2 2 

+TrC/TrC/ n - 1 VC/ G SU(7) (A5) 

TrU n = -TrU n - s + TrU 1 TrU n - 7 - — (Trf/) 6 Trf/"™ 6 + — (TrC/) 4 TrC/ 2 TrC/ n ~ 6 

720 V ; 48 V ; 

-^(Trf/) 2 (Tr[/ 2 ) 2 Trf/ n - 6 + ^(TrC/ 2 ) 3 TrC/ n - 6 - -^(TrC/) 3 TrC/ 3 TrC/ n - 6 

+^TiUTTU 2 TiU 3 TTU n - 6 - ^(TrC/ 3 ) 2 TrC/ n ~ 6 + i(Trf/) 2 Trf/ 4 Tr^ 6 

--TrU 2 TrU 4 TrU n - 6 - ^TrC/TrC/ 5 TrC/ n_6 + -TrC/ 6 TrC/ n_6 + — (TrC/flYC/™ -5 
8 5 6 120 v ; 

- — (TrC/) 3 TrC/ 2 TrC/ n_5 + ^TyU {TyU 2 ) 2 TrC/ n_5 + ^(TyU) 2 TyU 3 TyU 71 ' 5 

_lTrf/ 2 Trf/ 3 Trf/ n - 5 - iTrf/Tr^Trf/™- 5 + -TrC/ 5 TrCT~ 5 - — (TrC/) 4 TrC/ n - 4 
6 4 5 24 v ; 

+-(TrC/) 2 TrC/ 2 TrC/ n - 4 - -(TrU 2 ) 2 TrU n -' 1 - -TrUTrU 3 TrU n ~' 1 + -TrC/ 4 TYC/ n - 4 
4 V ; 8 V ; 3 4 

+i(TrCT) 3 TrC/ n ~ 3 - -TrC/TrC/ 2 TrC/ n_3 + -TrU 3 TrU n - 3 - -(TrC/flYC/™ -2 
6 2 3 2 

+-TYU 2 TrU n ~ 2 + TrC/TrC/ n_1 VC/ G SU(8) (A6) 
2 
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